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Abstract 

Underwater  acoustic  propagation  is  characterized  by  multipath  or  multimode  prop¬ 
agation.  Ray  theory  and  mode  theory  are  not  fully  adequate  for  modeling  physical 
reality.  Impulse  responses  can  be  more  accurately  calculated  using  Gaussian  beam 
theory.  Signal  processors  can  be  designed  to  take  advantage  of  the  channel  complex¬ 
ity  if  the  propagation  is  actually  known  so  that  detectability  is  increased.  The  pro¬ 
posed  technique,  channel  matched  filtering,  synthetically  back  propagates  the  wave 
front  to  a  hypothesized  source  location.  Accurate  passive  estimates  can  be  made 
without  knowledge  of  signal  characteristics.  GB  theory  can  easily  accommodate  a 
range-dependent  deep  water  environment. 


Introduction 


Underwater  acoustics  has  always  been  hampered  by  a  complicated  multipath  envi¬ 
ronment.  It  can  be  demonstrated  however,  that  the  complex  arrival  structure  can 
actually  be  used  to  the  signal  processor’s  advantage.  The  essential  element  is,  of 
course,  good  modeling  of  ocean  acoustic  properties.  The  study  of  techniques  which 
utilize  precise  knowledge  of  the  acoustic  propagation  is  the  proposed  thesis  problem. 

It  is  well  known  that  low  frequency  sound  (~  less  than  lKHz)  can  travel  several 
megameters  underwater  without  prohibative  attenuation  (Ref.  12).  The  transmis¬ 
sion  of  energy  over  such  great  distances  is  facilitated  by  two  facts.  The  first  is  that 
sound  is  attenuated  about  25dB  per  Mm  per  (kHz)2  at  low  frequencies.  This  means 
that  the  medium  does  not  absorb  too  much  of  the  sound  energy.  The  second  is 
that  the  sound  speed  profile  (fig.  1)  has  a  minimum  at  a  depth  of  approximately 
one  kilometer.  This  means  that  the  energy  is  ducted  as  in  a  waveguide  and  that 
boundary  losses  are  not  incurred  for  purely  refracted  rays. 

Another  well  documented  characteristic  of  the  ocean  acoustic  channel  is  multi- 
path  propagation  (Ref.  6).  A  typical  impulse  response  is  shown  in  figure  2.  More 
insight  into  the  physics  of  propagation  can  be  obtained  by  considering  the  timefront 
of  an  impulse.  Timefronts  are  the  locus  of  endpoints  of  rays  which  have  travelled  an 
equal  amount  of  time  from  the  source.  Figures  3-5  show  the  evolution  of  a  timefront. 
As  the  sound  propagates  further  and  further,  the  timefront  folds  accordion-like  gen¬ 
erating  new  sheets.  The  sheets  of  the  timefront  correspond  to  the  different  arrivals 
that  are  apparent  in  the  impulse  response. 

In  order  to  deal  with  the  multipath  propagation  problem  a  new  signal  processor 
is  proposed  which  synthetically  back-propagates  the  wavefront  to  a  hypothesized 
source  location.  This  technique  is  analagous  to  phase  conjugation  methods  popular 
in  optics.  It  is  similar  to  the  matched  field  processing  first  described  by  Bucker 
(Ref.  1)  and  more  recently  by  Clay  (Ref.  3).  Wax  and  Kailath  have  considered 
localization,  but  not  in  the  multipath  ocean  environment.  The  motivation  for  the 
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back-propagation  method  is  to  refocus  the  signal  energy,  which  is  dispersed  both 
in  time  and  space,  to  a  particular  source  location.  By  concentrating  the  energy, 
signal  detectability  should  increase.  By  searching  over  a  set  of  hypothesized  source 
locations,  a  source  position  estimate  can  be  made. 


Ocean  Acoustic  Modeling 


Good  modeling  of  ocean  acoustic  propagation  is  essential  to  good  signal  process¬ 
ing.  Efficiency  or  ease  of  calculation  is  just  as  important  as  model  accuracy  and 
consistency.  There  are  two  main  schools  of  thought  in  acoustic  modeling.  They 
are  known  as  ray  theory  and  mode  theory.  Neither  of  these  is  totally  acceptable 
for  our  purposes.  We  shall  consider  their  shortcomings  and  suggest  gaussian  beam 
theory  as  an  improvement  on  ray  theory.  In  order  to  simplify  this  discussion,  we 
shall  consider  only  two-dimensional  models  here. 

All  models  attempt  to  provide  a  solution  to  the  wave  equation. 

d2u  d2u  \_d2u 
dx 2  dz 2  c2  dt2 

In  this  formulation,  x  is  the  range  coordinate,  z  is  the  depth  coordinate,  c  =  c{x,z) 
denotes  the  sound  speed,  and  t  is  time.  The  quantity  u(x,  z,  t)  may  represent  various 
physical  quantities,  such  as  the  pressure  or  density. 

In  the  ray  theory  view  of  the  world,  one  assumes  a  solution  to  the  wave  equation 
of  the  form: 

u(x,z,t)  =  A(*,*)exp^‘-r<r-*>l 

Now  inserting  the  solution  into  the  wave  equation  and  taking  the  limit  as  u  goes  to 
infinity,  one  gets  the  eikonal  equation. 

(©■♦(©*-* 

The  locus  of  points  where  T  is  constant  defines  an  isotemporal  surface,  or  timefront. 

In  order  to  model  a  point  source  with  ray  theory  one  traces  a  number  of  rays  from 
the  source,  over  an  angular  distribution,  to  the  receiver.  The  amplitude  response  is 
then  inversely  proportional  to  the  spacing  of  adjacent  rays  at  the  receiver.  For  an 
accurate  estimate  of  arrival  time,  one  must  compute  the  travel  time  of  an  eigenray, 
the  ray  that  goes  through  the  receiver  location. 


Tracing  rays  on  a  digital  computer  is  a  fairly  easy  thing  to  do.  This  accounts 
for  the  popularity  of  ray  theory.  The  arrival  time  estimates  are  very  accurate  and 
it  is  easy  to  include  a  range  dependent  sound  speed  profile.  The  main  disadvantage 
with  ray  theory  is  that  the  amplitude  estimate  is  not  always  accurate.  The  theory 
predicts  infinite  amplitude  at  caustics,  which  are  where  neighboring  raypaths  cross. 
It  predicts  zero  amplitude  in  shadow  zones.  These  do  not  match  physical  reality. 
Another  problem  is  that  in  order  to  get  accurate  time  of  arrival  estimates  one  must 
trace  eigenrays.  Practically  speaking  this  means  that  many  rays  with  different 
initial  conditions  must  be  traced  until  a  ray  that  goes  through  (or  comes  close  to) 
the  receiver  location  is  found. 

Mode  theory  solutions  to  the  wave  equation  overcome  some  of  the  inadequacies 
of  ray  theory.  The  cost  is  complexity  of  calculation  and  a  very  small  collection  of 
boundary  geometries  that  admit  solutions.  The  solution  entails  the  calculation  of 
eigenfunctions  which  are  frequency  dependent.  A  necessary  assumption  is  that  the 
sound  speed  is  range  independent.  This  means  that  the  eigenfunctions  vary  with 
depth  only.  One  approach  to  the  solution  is  to  match  the  boundary  conditions 
at  the  surface  and  then  iterate  the  eigenvalue  until  the  boundary  conditions  at  the 
bottom  are  satisfied.  Given  the  correct  eigenfunctions,  the  amplitude  response  of  the 
channel  can  be  accurately  calculated  anywhere  in  the  channel.  So  for  a  broadband 
signal  the  eigenfunctions  must  be  calculated  across  the  spectrum.  With  all  these 
difficulties,  it  is  not  hard  to  understand  why  ray  theory  is  popular. 

There  is  an  alternative.  Gaussian  beam  theory  offers  the  best  of  both  theories. 
It  allows  easy  calculation  of  the  channel  impulse  response  and  accurate  amplitude 
estimates.  Also  it  can  use  a  range  dependent  sound  speed  profile  and  it  does  not 
require  eigenray  tracing  for  accurate  time  of  arrival  estimates. 

It  is  well  known  that  the  high  frequency  wave  field  propagates  mostly  along  rays. 
Gaussian  beam  theory  seeks  solutions  to  the  wave  equation  which  are  concentrated 
close  to  each  selected  ray  ft.  Define  an  orthogonal  coordinate  system  (s,  n)  along 
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the  ray.  The  coordinate  s  measures  the  arclength  along  the  ray  from  an  arbitrary 
reference  point,  n  represents  a  length  coordinate  in  the  direction  perpendicular  to 
ft  at  s.  Rewrite  the  wave  equation  in  the  new  coordinate  system.  The  parabolic 
approximation  to  the  wave  equation  is  used  to  obtain  solutions  close  to  the  ray.  A 
complete  derivation  is  given  by  Cerveny  (Ref.  2).  The  solution  which  has  the  form 

/  t  r^)l1/2  f.  f3  ds  \  iu  p(s)  2 

= 4  ta)]  exp  b“  ('  -  L  «r>) +  Tiw”  ] 

can  then  be  found. 

The  central  axis  of  the  beam  obeys  the  standard  ray  trace  equations.  For  the  two 
dimensional  case  a  system  of  four  first  order  linear  differential  equations  describes 
the  raypaths. 

dt 


dt 


A  and  B  are  dummy  variables.  The  initial  conditions  are  related  to  the  source 
location  and  the  initial  angle  that  the  ray  has  when  it  leaves  the  source.  The  values 
of  p  and  q  along  the  ray  can  then  be  found  as  the  solution  to  two  additional  first 
order  linear  differential  equations  subject  to  some  initial  conditions. 


2“  =  c{s)p(s) 

Where  c„n  denotes  the  second  partial  derivative  of  the  sound  speed  with  respect  to 
the  normal  coordinate.  Both  p  and  q  are  complex  quantities  and  do  not  depend  on 
frequency.  The  solution  can  be  rewritten  as 

“<s' 1 0  =  *  [jwl  '  exp  -  r<s>l +  fcK{,) +  z£j 
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where 


.»  «■>  ‘.l 


*’<•)- <«*•({«) 
L(a) = [l/m  (M)} 
™=££) 


From  the  physical  point  of  view,  K  denotes  the  curvature  of  the  phase  front  of  the 
beam,  and  £  is  a  frequency  dependent  effective  half-width  of  the  beam. 

Having  solved  the  parabolic  wave  equation  near  a  single  beam,  if  we  want  to 


model  a  source  we  must  fit  a  set  of  Gaussian  beams  to  the  initial  wave  field.  The 


wave  field  u(M)  at  the  point  M  can  be  expressed  in  terms  of  Gaussian  beams  as 


u(M)  =  £$(<6)u*(s,n) 


where  <p  is  the  index  connected  with  ray  ft,  and  the  arbitrary  function  $(<£)  is 
specified  by  the  initial  source  field. 

Once  the  wave  field  due  to  a  source  at  point  S  and  a  receiver  at  any  point  7 Z 
in  the  ocean  acoustic  channel  is  known,  the  impulse  response  from  S  to  72  can  be 


written. 


M0  =  £>M(<-n) 


The  advantage  of  Gaussian  beam  theory  is  that  for  only  slightly  more  computational 
effort  than  for  conventional  ray  theory,  accurate  values  for  the  A  ’s  can  be  calculated. 
It  should  be  pointed  out  that  although  p  and  q  are  independent  of  frequency,  A,  and 
Ti  are  theoretically  weakly  dependent  on  frequency.  Our  best  information  is  that  in 
the  real  ocean  Ai  and  ry  do  not  change  appreciably  over  the  range  of  frequencies 
(5-500  Hz)  that  we  are  interested  in  (Ref.  5).  The  question  of  sensitivity  of  the 
impulse  response  to  frequency  needs  to  be  studied  further. 

All  the  models  considered  so  far  are  deterministic.  Ocean  variability  is  a  very 
complex  phenomena  but  its  main  component  is  thought  to  be  internal  waves.  Ex¬ 
perimentally  it  has  been  shown  that  it  is  the  amplitude  of  the  arrivals  that  varies 


the  most.  Amplitude  stability  is  on  the  order  of  ten  minutes.  Time  of  arrival  is 
actually  quite  stable,  variations  less  than  milliseconds  over  periods  of  days  (Ref.  6). 
E%en  crude  simulations  of  internal  wave  effects  confirm  this  behavior.  Investigation 
of  amplitude  variance  vis-a-vis  travel  time  variance  would  be  reasonable  and  inter¬ 
esting.  It  has  been  hypothesized  that  d2c/dz 2  affects  the  former  strongly  and  the 
latter  weakly. 
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Signal  Processing 


Taking  a  systems  approach,  we  want  to  apply  statistical  communications  theory  and 
signal  detection  and  estimation  theory  to  the  ocean  acoustic  situation.  It  is  crucial 
to  realize  that  this  is  a  multiple  source  and  multiple  sensor  environment  and  that 
sources  and  sensors  are  coupled  by  more  than  one  path. 

The  question  we  wish  to  answer  is  this:  If  we  knew  c(x,  y ,  z,  t),  the  speed  of  sound 
everywhere,  what  good  would  it  do?  Would  we  be  able  to  use  the  knowledge  of  the 
acoustic  transfer  function  to  improve  communications,  detection,  or  localization? 

We  choose  to  focus  on  signal  detectability  first  and  we  do  not  wish  to  make 
any  assumptions  about  the  nature  of  the  source  signals.  Before  proposing  any  new 
detection  algorithms  one  has  to  ask  why  is  the  standard  theory  of  signal  detectabil¬ 
ity  inadequate.  The  likelihood  ratio  theory  of  signal  detectability  is  essentially  a 
hypothesis  testing  theory  (Ref.  1,  13).  The  most  general  problem  might  be  set  up 
by  considering  a  distribution  of  source  locations  in  three  dimensional  space.  Given 
a  specific  realization  72  of  Af  source  locations,  one  specifies  a  composite  distribution 
for  each  source  waveform.  There  is  a  conditional  likelihood  ratio  l{ observation  |  72) 
with  respect  to  some  measure.  The  posterior  probability  of  a  particular  source 
distribution  would  entail  an  expectation  over  all  72. 

Detection  theory  does  give  some  guidelines.  It  suggests  highly  parallel  processing 
of  the  observation,  one  branch  for  each  hypothesis  (perhaps  each  acoustic  source). 
This  in  turn  suggests  that  linear  preprocessing  be  considered  just  prior  to  the  in¬ 
evitable  non-linear  processing  that  precedes  expectations  over  signal  parameters  and 
final  decisions. 

Detection  theory  uses  a  forward  propagating  approach:  it  hypothesizes  the  causes 
(source  positions  and  waveforms)  and  calculates  how  likely  each  would  be  to  have 
caused  the  observed  waveform.  Given  the  muliplicity  of  source  locations,  sensor 
positions,  and  paths,  this  appears  to  be  a  monumental  task.  Therefore  a  direct 
global  attack  such  as  this  will  not  be  undertaken  and  alternatives  should  be  sought 


that  do  not  need  to  hypothesize  the  waveforms  of  sources. 

As  shown  in  the  last  section,  the  ocean  acoustic  channel  can  be  modeled  as  linear 
time-invariant  filter  of  the  form 

P 

hU(T)  =  J2ap6(t  -  rp) 

p= 1 

where  U  represents  the  position  of  the  hydrophone  relative  to  the  source.  We  will 
use  this  deterministic  model  and  ignore  any  variability  of  the  ocean  since  it  is  our 
objective  to  develop  new  signal  processing  strategies  based  on  a  known  acoustic 
transfer  function. 

The  new  processing  strategy  we  propose  to  study  is  the  reverse  propagation 
of  signal  waveform  through  a  linear  filter  which  approximates  the  ocean  channel. 
Suppose  we  wish  to  find  the  source  waveform  from  the  waveform  observed  at  the 
hydrophone.  We  will  assume  that  the  observation,  r(<)  =  s(t)  *  hu(t),  is  noiseless 
for  the  time  being.  The  reverse  propagating  filter  has  the  form 

Q 

bv(T)  =  '£B,6(T  +  r<l) 

where  V  is  the  hypothesized  position  of  the  source  relative  to  the  hydrophone,  and 
the  Bq's  and  r,’s  are  related  to  the  forward  coefficients  and  time  delays  respectively. 
After  the  filter,  the  output  is 

a  P  Q 

w(t)  =  r(t)  *  6v(t)  =  "  Tv  +  Ti) 

p=l  9=1 

Now  if  the  reverse  propagating  filter  has  correctly  hypothesized  the  source  position, 
w(t)  becomes 

w(t)  =  £  ApBps(t)  +  £;  £;  ApBqs(t  -  TP  +  r, ) 

p=  i  p=1  ?#p 

The  first  term  has  the  property  we  want  of  reconstructing  the  original  unknown 
source  waveform.  The  second  term  is  a  self-interference  term. 


Let  us  now  examine  how  we  might  choose  the  coefficients,  Bq,  of  the  reverse 
propagating  filter.  Three  choices  come  to  mind. 

Bq  =  A\ 

where  z  —  1.0.  or  —1.  These  correspond  to  a  channel  matched ,  a  channel  phase ,  and 
a  channel  inverse  filter  respectively.  As  z  varies  between  1  and  —1,  we  might  expect 
that  we  could  trade  signal-to-noise(or  interference)  ratio  for  sharpness  of  response. 

If  one  knew  the  signal  and  noise  power  spectral  density  at  the  hydrophone,  the 
matched  filter  for  the  received  waveform  would  be 

\N{u)\ ^ 

The  channel  matched  filter  is  nothing  but  the  channel  half  of  the  matched  filter. 

The  channel  phase  and  channel  inverse  filters  have  a  less  clear  relationship  to  the 
standard  phase  only  and  inverse  filters. 

Next,  let  us  assume  that  we  have  an  array  of  hydrophones.  Each  hydrophone 
reception  is  back-propagated  to  a  hypothesized  source  location  and  then  the  filter 
outputs  are  summed  to  get 

a(r;  tf,  V)  =  £  £  £  APAV)BqAV)6(r  -  rp,r(f/)  +  r,,r( V)) 

p—  1  1  T—  1 

where  r  is  the  hydrophone  index  and  R  is  the  number  of  hydrophones.  The  function 
a(r;  U,  V)  is  a  psuedo-ambiguity  function.  It  can  be  expected  to  peak  when  U  =  V, 
but  the  volume  enclosed  by  it  is  not  constant  as  V  varies. 

Let  us  assume  that  we  are  using  channel  matched  filtering.  There  are  two  cases 
which  can  occur.  First  suppose  that  the  back-propagating  filter  exactly  compensates 
for  the  channel.  Then  we  have 

PR  P  Q  R 

a(t;  U,  V)  =  X;2>P,r(C0A>,r(W0+L  EEApAU)AqAUMt-rvAU)+TqAU)) 

P=l r=1  p=l 9#pr=l 

The  first  term  reconstructs  the  original  signal  and  the  second  is  a  self-interference 
term.  The  hope  is  that  as  R  and  P  increase,  the  gain  applied  to  the  reconstruction 


increases  faster  than  that  of  the  interfering  term.  This  depends  upon  the  unknown 
power  spectral  density  of  the  source  signal.  Now  suppose  that  the  back -propagating 
filter  incorrectly  hypothesizes  the  source  position.  Then  we  have 

P  Q  R 

a(t;U,V)  =  '£EYi^AU)A1AV)s(t-rJ>AU)  +  TqAV)) 

p=  1  ij  =  l  r  =  1 

In  this  case  we  have  only  the  self-interference  term.  Adding  many  signal  waveforms 
in  this  incoherent  fashion,  we  hope  that  the  power  in  a(t;U,V)  is  less  than  in 
a(t\U,U )  where  the  waveforms  are  added  coherently. 

It  may  be  useful  to  perform  a  normalization 

a(t 

This  assures  that  |a|  <  1  by  an  application  of  the  Schwarz  inequality.  Further 
processing  of  the  pseudo-ambiguity  function  would  probably  resemble  standard  de¬ 
tection  theory.  For  example  the  integral 

\A(U,V)\2±  I Ta\t;U,V)dt 
Jo 

could  be  formed  and  an  energy  detector  would  then  be  realized. 

The  ability  of  an  array  of  back-propagating  filters  to  localize  a  source  depends 
on  the  complexity  of  the  arrival  structure.  Generally  speaking  the  more  complicated 
the  structure  the  easier  it  becomes  to  localize  a  source.  This  happens  because  of 
the  fact  that  more  structure  is  equivalent  to  more  information  about  the  channel. 
Thus  we  have  used  the  multipath  channel  characteristic  to  our  advantage. 

The  total  array  signal  gain  is  the  ratio  of  the  power  in  the  reconstructed  signal 
term  to  the  power  of  the  signal  in  the  largest  path.  The  power  in  the  reconstructed 
signal  term  is  proportional  to  PR,  the  number  of  hydrophones  R,  and  the  number 
of  paths  P.  One  cannot  say  what  the  reconstruced  signal  power  will  be  in  gen¬ 
eral  without  specifing  where  the  source  is  because  the  answer  depends  on  the  path 
coefficients  «4p,r-  The  reason  for  choosing  the  power  in  the  largest  path  for  the  de¬ 
nominator  instead  of  the  power  in  all  paths  is  this:  In  conventional  signal  processing 


;U,V)  = 


a(t;U,V) 

~  |B,.r(VOIJ  “  |r(t  +  r,,,(r))|! 
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one  cannot  recombine  paths  and  therefore  one  views  the  other  paths  as  contributing 
interference. 

The  noise  performance  of  our  system  is  a  hard  thing  to  estimate.  This  is  be¬ 
cause  the  noise  characteristics  of  the  ocean  are  not  known  well.  So  much  of  what 
is  normally  refered  to  as  noise  is  actually  interference  from  other  sources  such  as 
shipping  traffic,  active  sonars,  and  drilling  rigs.  If  it  is  possible  to  map  all  the  major 
sound  sources  in  the  ocean  then  much  of  this  interference  could  be  eliminated  and 
the  ultimate  ocean  noise  limit  could  be  found. 


Results 


Results  to  date  have  been  very  encouraging,  both  in  terms  of  modeling  and  in  terms 
of  signal  processing. 

Impulse  responses  have  been  calculated  on  a  computer  using  Gaussian  beam 
theory.  This  is  new  work  and  we  do  not  believe  anyone  else  has  calculated  point- 
to-point  transfer  functions  using  this  theory.  Although  it  has  not  been  verified  by 
comparison  with  wave  theory,  the  responses  exhibit  the  correct  behavior  at  caustics 
and  in  shadow  zones.  An  example  impulse  response  between  two  points  on  the 
SOFAR  axis  is  given  in  figure  2.  This  was  calculated  using  a  Munk  canonical  sound 
speed  profile  and  the  two  points  are  500  kilometers  apart. 

The  pseudoambiguity  function  has  been  calculated  in  a  number  of  test  cases  and 
it  exhibits  the  correct  response.  The  response  peaks  at  the  true  source  location  and 
decreases  as  the  hypothesized  source  position  moves  away  from  the  true  position. 
In  figure  6,  the  maximum  over  time  of  the  unnormalized  matched  pseudoambiguity 
function  is  plotted  versus  range. 

max a(r;  U,V)  =  maxf^f^ PAU)A,AV)6(r  -  rPAU)  +  rq,r(V)) 

p=l  fsl  r=l 

The  true  source  location  is  at  250  kilometers,  where  the  peak  is.  The  search  depth 
was  held  fixed  at  the  true  search  depth.  Channel  matched  filtering  works. 

For  the  next  several  cases  we  will  consider  a  slightly  different  function.  Suppose 
the  source  waveform  is  known  to  be  a  sinusoid.  Then  the  pseudoambiguity  function 
becomes 

P  Q  R 

a(t;  U,  V)  =  £  £  5>p,r((/)S7,r(ncos[u;(t  -  rPAU)  +  rt,,(V))] 

p=l  7=1  r=  1 

We  will  evaluate  this  function  at  t  =  0  for  a  variety  of  cases. 

Consider  first  the  phase  only  filter.  The  response  is  plotted  versus  search  range 
for  a  source  at  500  Hz  in  figures  7-9.  The  function  becomes  better  behaved  as  the 
true  range  is  increased  and  the  number  of  hydrophones  in  the  array  is  increased. 


Now  we  will  consider  the  matched  filter  case.  In  figure  10,  the  source  is  located 
100  kilometers  from  the  array.  Both  source  and  array  are  at  a  depth  of  1200  meters. 
The  response  function  becomes  much  smoother  than  in  the  phase  only  filter  case. 
The  sidelobe  structure  also  becomes  much  more  defined.  The  magnitude  of  the  side- 
lobes  is  a  cause  for  worry  though.  It  is  postulated  that  this  behavior  occurs  because 
of  convergence  zone  focusing.  In  figure  11,  the  number  of  phones  is  increased  to  500. 
This  has  the  effect  of  reducing  the  size  of  sidelobes  compared  to  the  unnormalized 
main  response.  The  true  source  location  is  moved  to  250  kilometers  in  figure  12. 
This  also  has  the  effect  of  reducing  the  sidelobe  response  even  though  in  this  case 
only  50  phones  are  used  in  the  array.  This  is  to  be  expected  because  of  the  greater 
complexity  of  the  wavefront  as  range  increases. 

Let  us  now  examine  the  depth  resolution  of  the  processing.  In  figures  13  and  14, 
the  source  depth  is  400  meters  and  the  array  depth  is  1200  meters.  The  search  depth 
is  400  meters  in  figure  13  and  1200  meters  in  figure  14.  The  response  peaks  at  the 
correct  location  in  figure  13.  When  the  search  depth  is  changed  the  absolute  value 
of  the  unnormalized  response  is  much  lower.  The  peak  at  about  290  kilometers  in 
both  plots  was  unexpected  and  is  so  far  unexplained. 

Finally,  the  last  plot,  in  figure  15,  shows  the  effect  of  a  frequency  change  in  the 
pseudoambiguity  function.  The  parameters  are  the  same  as  in  figure  13  except  the 


Future  Work 


As  seen  in  the  last  section,  channel  matched  filtering  is  a  viable  technique.  Many 
basic  questions  about  the  performance  of  back-propagating  processing  need  to  be 
answered  however.  The  performance  measures  which  we  are  interested  in,  include 
signal  detectability,  localization  resolution,  array  gain,  and  noise  or  interference 
rejection.  The  two  most  critical  issues  are  probably  array  geometry  and  location 
aliasing. 

Array  geometry  is  somewhat  limited  by  physical  constraints.  We  wish  to  examine 
the  theoretical  performance  of  different  configurations  without  regard  to  hardware  or 
computational  costs  however.  Different  configurations  could  include:  a  small  densely 
instrumented  sensor  array  in  the  horizontal;  a  convergence  zone  sized  sparsely  instru¬ 
mented  sensor  array  in  the  horizontal;  a  top-to-bottom  stiff  densely  instrumented 
vertical  array  in  conjunction  with  a  conventional  sparsely  instrumented  horizontal 
array. 

Location  aliasing  is  manifested  in  the  ambiguity  of  which  convergence  zone  a 
source  may  be  in.  This  question  is  probably  intimately  related  to  the  particular  array 
geometry.  Performance  in  this  respect  could  possibly  be  improved  by  implementing 
a  scheme  similar  to  null  steering;  that  is,  by  trying  to  reduce  the  sidelobe  response 
where  it  is  magnified  by  convergence  zone  focusing.  A  full  scale,  three  dimensional 
examination  of  the  pseudoambiguity  function  would  be  informative.  The  increase 
in  dimensionality  might  improve  localization  performance  since  would  increase  the 
complexity  of  the  wavefronts  considered. 

Another  important  issue  is  the  signal  characteristics.  As  signal  bandwidth  is 
increased  the  performance  of  the  array  should  improve.  This  happens  because  more 
bandwidth  would  again  increase  the  waveform  complexity. 

One  question  is  what  is  the  ultimate  performance  limit  of  matched  channel 
processing.  A  useful  approach  may  be  to  consider  how  well  we  could  image  a  source 
if  we  know  the  wavefront  exactly  on  some  continuous  smooth  surface.  Perhaps  we 
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could  learn  what  the  best  surface  sheet  geometry  would  be.  Also  it  might  be  possible 
to  relate  resolution  to  aperture  size. 

If  two  or  more  sources  are  present,  standard  techniques  in  signal  processing  may 
be  applied  to  improve  performance.  In  order  to  seperate  sources  high  resolution 
spectral  estimation  techniques  may  be  used  (Ref.  9).  Also  recursive  methods  (Ref. 
10)  could  be  employed  to  cancel  energy  from  known  or  previously  identified  inter¬ 
feres. 

An  improvement  may  be  made  to  the  channel  matched  filter  processing  by  at¬ 
tempting  to  resolve  the  paths.  The  Gaussian  beam  model  specifies  the  angle  of 
arrival  of  different  paths.  This  information  is  ignored  currently  and  the  waveform 
is  back-propagated  in  all  directions.  If  one  is  able  to  resolve  the  paths  (which  are 
distibuted  over  about  30  degrees)  then  the  back-propagation  could  be  done  only  in 
the  direction  of  the  hypothesized  source.  This  technique  could  vastly  improve  the 
self-inter  -rence  performance. 

Finding  the  location  of  a  source  by  channel  matched  filtering  can  be  looked  upon 
as  the  solution  of  an  inverse  problem  by  forward  modeling.  In  the  usual  situation 
we  have 

R(f)  =  S(f)H(f)  +  N(f) 


where  S(f)  is  the  source  spectrum,  H(f)  is  the  channel  transfer  function,  and  R(f)  is 
the  observed  spectrum.  We  wish  to  form  an  estimate  of  some  channel  parameters 
such  as  the  range,  depth,  and  bearing  angle.  This  task  may  be  performed  in  two 
ways.  The  usual  method  is  to  directly  perform  an  inversion  operation  on  R(f) 


*(/)* 


*(/) 

S(f) 


to  obtain  the  channel  parameters.  The  problem  with  this  method  is  that  the  in¬ 
version  operation  is  often  very  sensitive.  That  is  to  say  small  changes  in  R(f)  will 
produce  large  changes  in  the  estimated  channel  parameters.  This  is  directly  solving 
the  inverse  problem. 
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There  is  another  method  of  solving  for  the  channel  parameters.  In  this  method 
a  model  for  H(f)  is  postulated  that  uses  the  desired  estimates  of  channel  parameters 
as  inputs.  A  whole  class  of  transfer  functions  can  then  be  generated.  Note  that  this 
class  is  a  subset  of  the  class  of  transfer  functions  estimated  by  the  direct  procedure. 
The  best  estimate  of  the  channel  parameters  is  then  found  by  finding  the  H(f) 
which  maximizes 

R(f)S’(f)H’(f) 

l*V(/)|2 

This  is  the  structure  that  is  used  in  the  proposed  thesis  problem  and  is  known  as 
solving  the  inverse  problem  by  forward  modeling.  Formalizing  this  algorithm  would 
be  an  interesting  and  useful  exercise.  Some  work  has  been  done  in  this  area  by 
Jaynes  (Ref.  8). 

Another  consideration  is  the  effect  of  Doppler  shifts  on  the  processing.  When  an 
acoustic  source  or  receiver  is  moving  with  respect  to  the  medium  a  frequency  shift 
in  the  received  signal  is  observed.  Note  that  this  shift  is  different  for  all  paths.  This 
is  because  all  paths  have  different  takeoff  and  arrival  angles.  The  amount  of  doppler 
shift  is  then  proportional  to  the  dot  product  of  the  source  (or  receiver)  velocity 
vector  and  the  path  takeoff  (or  arrival)  vector. 

Finally,  the  application  of  channel  matched  filtering  to  real  data  must  be  con¬ 
sidered.  How  this  might  be  done  is  not  dear.  There  is  much  data  available  on  long 
range  acoustic  transfer  functions.  The  data  was  generated  by  the  Ocean  Acoustic 
Tomagraphy  Group  of  which  CSPL  is  a  member.  This  data  could  therefore  be  easily 
accessed. 
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Figures 

Figure  i  is  Munk's  canonical  sound  speed  profile. 

c{z)  =  c0[l  +  £(expr>  +77  —  1 )] 
r>  =  2(z  -  zo)/0 

This  is  a  theoretical  profile  developed  from  first  order  physical  reasoning,  and  rep¬ 
resents  an  oversimplified  model  of  real  sound  speeds. 

Figure  2  was  calculated  for  a  range  of  500  kilometers,  a  source  depth  of  1200 
meters,  and  a  receiver  depth  of  1200  meters.  The  times  marked  correspond  to  the 
time  of  propagation  at  the  edges  of  the  graph.  The  height  of  the  impulses  are 
proportional  to  the  amplitude  of  the  response. 

Figures  3-5  show  stages  of  the  timefront  for  a  source  depth  of  800  meters  and 
a  Munk  sound  speed  profile.  The  horizontal  axis  is  range,  but  magnified  by  1.5 
compared  to  vertical.  Purely  refracted  rays  with  source  angles  between  ±15  degrees 
are  shown. 

Figure  3:  After  1  second  the  wavefront  is  nearly  spherical. 

Figure  4:  After  68  seconds,  two  sets  of  wavefronts  have  developed.  They  essen¬ 
tially  correspond  to  initially  upward  going  and  downward  going  rays. 

Figure  5:  After  337  seconds,  the  timefront  is  the  interleave  of  two  accordions, 
each  with  8  sheets.  A  receiver  at  a  depth  of  1200  meters  (on  the  axis)  will  have  the 
time  arrival  structure  shown  in  figure  2,  with  groups  of  4  arrivals. 

Figure  6-15  assume  a  horizontal  line  array  with  variable  phone  spacing  and 
that  the  source  is  placed  in  an  endfire  position  relative  to  the  array.  All  ambiguity 
functions  were  calculated  assuming  a  range  independent  Munk  profile. 

Figure  6:  Maximum  of  ambiguity  function  with  impulses  sorted  into  1  millisec¬ 
ond  time  bins.  True  source  range  250,000  meters,  true  source  depth  1200  meters, 
99  phones  spaced  100  meters  apart  at  1200  meters  depth,  matched  processing  at  a 
search  depth  of  1200  meters. 


Figure  7:  Ambiguity  function  with  sinusoidal  source  at  500  hertz,  true  source 
range  100,000  meters,  true  source  depth  1200  meters,  100  phones  spaced  100  meters 
apart  at  1200  meters  depth,  phase  only  processing  at  a  search  depth  of  1200  meters. 

Figure  8:  Ambiguity  function  with  sinusoidal  source  at  500  hertz,  true  source 
range  250,000  meters,  true  source  depth  1200  meters,  500  phones  spaced  100  meters 
apart  at  1200  meters  depth,  phase  only  processing  at  a  search  depth  of  1200  meters. 

Figure  9:  Ambiguity  function  with  sinusoidal  source  at  500  hertz,  true  source 
range  500,000  meters,  true  source  depth  1200  meters,  200  phones  spaced  100  meters 
anart  at  1200  meters  depth,  phase  only  processing  at  a  search  depth  of  1200  meters. 

Figure  10:  Ambiguity  function  with  sinusoidal  source  at  500  hertz,  true  source 
range  100,000  meters,  true  source  depth  1200  meters,  100  phones  spaced  100  meters 
apart  at  1200  meters  depth,  matched  processing  at  a  search  depth  of  1200  meters. 

Figure  11:  Ambiguity  function  with  sinusoidal  source  at  500  hertz,  true  source 
range  100,000  meters,  true  source  depth  1200  meters,  500  phones  spaced  100  meters 
apart  at  1200  meters  depth,  matched  processing  at  a  search  depth  of  1200  meters. 

Figure  12:  Ambiguity  function  with  sinusoidal  source  at  500  hertz,  true  source 
range  250,000  meters,  true  source  depth  1200  meters,  50  phones  spaced  100  meters 
apart  at  1200  meters  depth,  matched  processing  at  a  search  depth  of  1200  meters. 

Figure  13:  Ambiguity  function  with  sinusoidal  source  at  500  hertz,  true  source 
range  250,000  meters,  true  source  depth  400  meters,  100  phones  spaced  100  meters 
apart  at  1200  meters  depth,  matched  processing  at  a  search  depth  of  400  meters. 

Figure  14:  Ambiguity  function  with  sinusoidal  source  at  500  hertz,  true  source 
range  250,000  meters,  true  source  depth  400  meters,  100  phones  spaced  100  meters 
apart  at  1200  meters  depth,  matched  processing  at  a  search  depth  of  1200  meters. 

Figure  15:  Ambiguity  function  with  sinusoidal  source  at  100  hertz,  true  source 
range  250,000  meters,  true  source  depth  400  meters,  100  phones  spaced  100  meters 
apart  at  1200  meters  depth,  matched  processing  at  a  search  depth  of  1200  meters. 
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